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SUMMARY 

A series of test flights was conducted by the U. S. Navy over a 3- 
year period to evaluate the effects of icing on the operation of the 
ZPG-2 airship. In supercooled clouds, ice formed only on the forward 
edges of small protuberances and wires and presented no serious hazard 
to operation. Ice accretions of the glaze type which occurred in condi- 
tions described as freezing drizzle adversely affected various components 
to a somewhat greater extent. The results indicated a need for protec- 
tion of certain components such as antennas, propellers, and certain 
parts of the control system. 

The tests showed that icing of the large surface of the envelope 
occurred only in freezing rain or drizzle. Because of the infrequent 
occurrence of these conditions, the potential maximum severity could not 
be estimated from the test results. The increases in heaviness caused 
by icing in freezing rain and drizzle were substantial, but well within 
the operational capabilities of the airship. 

In order to estimate the potential operational significance of icing 
in freezing rain, theoretical calculations were used to estimate: (l) 

the rate of icing as a function of temperature and rainfall intensity, 

(2) the climatological probability of occurrence of various combinations 
of these variables, and (3) the significance of the warming influence of 
the ocean in alleviating freezing-rain conditions. The results of these 
calc u l a tions suggest that, although very heavy icing rates are possible 
-in combinations of low temperature and high rainfall rate, the occurrence 
of such conditions is very infrequent in coastal areas and virtually im- 
possible 200 or 300 miles offshore. 


INTRODUCTION 

The nonrigid airship ("blimp") has been adapted to carry aircraft 
detection equipment for use in the air defense system. The aircraft- 
early-warning mission requires that picket aircraft endure all types of 
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weather conditions while maintaining a continuous watch from a particular 
area. In support of this mission a project was established by the Office 
of Naval Research to evaluate the all-weather capabilities of the ZPG-2 
airship (975,000 cu ft envelope) (fig. l). The flight testing was assigned 
to the Naval Air Development unit at the United States Naval Air Station 
in South Weymouth, Massachusetts. The MCA Lewis laboratory was requested 
by the Office of Naval Research to assist the Navy by supplying instru- 
mentation and technical assistance in analyzing and studying the icing 
problem as applied to airships. 

Weather conditions considered hazardous to airborne operations are 
those conducive to ice and snow accumulations on various parts of the 
airship and those affecting the operational capability of the airship, 
such as winds and turbulence. The flight tests and study reported herein 
were conducted only for the purpose of investigating the airworthiness of 
the airship in various types of icing conditions. 

Previous airship operations in patrol and convoy duty were not seri- 
ously confronted with the icing problem because these missions were flown 
at altitudes (50 to 500 ft) below the levels of most icing clouds. 

Icing and snow problems of most concern for early airships were mostly 
confined to difficulties encountered during outside mooring on the mast. 
The added weight of ice and snow accumulations on the top of the envelope 
would buckle the landing gears and in cross winds would cause the airship 
to roll to one side and damage outriggers and fins. The in-flight ice 
and snow problem was investigated to a limited extent in 1945 when the 
MCA assisted the Navy in a preliminary study of ice-protection require- 
ments for the K-type airship (435,000 cu ft envelope). The results are 

reported in reference 1 for a 2-|-hour flight in conditions described as 

being a mixture of wet snow and freezing rain and which produced no ad- 
verse effects on the operation of the airship. The icing and snow prob- 
lem has become relevant to the present ZPG-2 airship with its increased 
altitude range in the aircraft -early-warning mission which places oper- 
ations more frequently at levels of icing clouds. 

The flight tests reported herein were conducted over a period of 
three icing seasons along the coast of the northeastern United States. 

A completely unprotected ZPG-2 airship was flown in icing conditions 
during the first two winters. A second similar airship incorporating 
limited icing protection based on experience gained from the earlier 
flights was used during the final season. Test objectives were to 
determine the components of the airship which collect ice and the result- 
ing effects on operation and performance . Flights were also planned to 
evaluate the ice and snow problem over water. The first part of this 
report describes several flights in a variety of ice and snow conditions 
and the resulting effects on various components as reported by Naval Air 
Development unit personnel. 
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A preliminary study of the airship icing problem indicated that 
accumulations of ice or snow on the large surface areas along the top 
of the envelope and fins could produce the serious problem of weight 
greater than the lifting capacity of the airship. This problem would 
exist only in the more unusual conditions of freezing drizzle or rain 
where the drop sizes are of sufficient size to strike the envelope sur- 
face. Droplets in the smaller size range of supercooled clouds (<100 
microns) would not impinge on the envelope in sufficient amounts to be 
of concern. The low airspeed and gradual curvature of the large air- 
ship envelope allow the small light-weight droplets to be easily de- 
flected by the airstream around the envelope. Atmospheric conditions 
encountered during the flight tests did not provide sufficient informa- 
tion to evaluate this primary problem of ice loads from freezing rain 
adequately. It was considered advisable, therefore, to calculate rates 
of ice formation on the airship envelope based on theoretically derived 
droplet trajectories and heat- transfer relations combined with rainfall 
rates and temperatures typical of freezing-rain conditions. The opera- 
tional significance of this problem during offshore missions was also 
considered desirable. The second part of this report presents the re- 
sults of these calculations, including an analysis of the probability 
of freezing rain for offshore areas. 


DESCRIPTION C3? TEST AIRSHIP AND INSTRUMENTATION 

The ZPG-2 airship used in the tests (fig. l) has a fabric envelope 
filled with approximately 975,000 cubic feet of helium. The over-all 
length is 342 feet and the maximum diameter is 75 feet. Two 18-foot 
propellers are mounted on outriggers projecting from each side of the 
car. Power is transmitted through shafting from the engines mounted in- 
side the car. The stabilizers and control surfaces (ruddevators ) are 
mounted at 45° angles and are supported by brace wires attached to the 
envelope surface. Exposed control cables run from the car to the rudde- 
vators. The airship is normally operated in a heavy condition (i.e., 
heavier than the air displaced by the envelope) and thus requires same 
dynamic lift for flight. 

The test instrumentation consisted of equipment to measure and ob- 
serve ice accretions on various airship components and to evaluate the 
effects of the ice on the operation of the airship. Rate- of -accretion 
meters indicated the magnitude of ice formations on the forward edges 
of small exposed components. These meters recorded the rate at which 
ipe formed an the leading edge of a thin rotating disk and were used in 
flight to detect the presence of icing and the duration of the icing 
conditions (see ref. 2 for a description of these meters). Ice deposited 
per unit area an the top surfaces of the envelope in freezing rain or 
drizzle was measured at three locations (nose, center, and tail sections) 
by use of a devipe developed specifically for this purpose by the Clevite 
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Brush Company. This instrument utilized the change in resonant frequency 
caused by ice accumulating on a vibrating sensing element mounted flush 
with the surface. The frequency change of piezo-electric driving stacks 
contained in the sensing elements was measured electronically. Onl y pre- 
liminary development of this instrument was accomplished during the pro- 
gram. The rotating -multicylinder technique was employed during part of 
the program to measure droplet sizes and liquid-water content of the icing 
clouds. Ice formations on areas not visible from the car were observed 
with movie cameras mounted at several points throughout the airship and 
with a television camera installed on the top to view the aft section. 

Flight conditions were either recorded on film or observed by the 
crew. Outside air temperature, altitude, and airspeed were continuously 
recorded. The flying qualities of the airship under test conditions were 
judged mainly by the pilot. For most flights the handling characteris- 
tics and increased load caused by the ice accretions were estimated by 
relating angle of attack, airspeed, and power settings to those required 
for normal ice-free operation. 


DISCUSSION OF ICING FLIGHTS 

Icing conditions of sufficient intensity and duration to provide 
pertinent results were encountered on nine airship flights along the 
coast of the northeastern United States during three icing seasons. 

Table I summarizes the data from these flights using the records from 
the disk-type icing-rate meters to evaluate the icing conditions. Rime 
ice was encountered on six flights and freezing rain or drizzle was 
present during three icing encounters. The rime ice, resulting from 
small supercooled cloud droplets, formed only on exposed objects and 
protuberances and not on the large envelope surface. The freezing- 
drizzle conditions, where the droplet sizes were much larger, produced 
a film of ice on the top envelope surfaces as well as glaze -ice forma- 
tions on exposed components. Compared with the rime- ice encounters, the 
freezing drizzle produced somewhat more deleterious effects on the oper- 
ation of the airship. The many factors, including droplet size, which 
influence the ice collection characteristics of various sized objects 
are discussed in a later section. 


Flights in Freezing Rain and Drizzle 

The largest total amount of ice measured during the program was 
accumulated on flight number 3. Glaze- ice formations were produced 
from a condition described as freezing drizzle. Although droplet sizes 
were not measured, the flight reports and ground observations indicated 
the airship was flown in clouds in an area of light precipitation. 
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The increasing ice accumulation with time as measured by the icing- 
rate meters is shown in figure 2 for a chronology of events reported 
during the flight . Icing was encountered during climb shortly after 
takeoff from South Weymouth. The greatest icing rate appeared to be at 
the 3000-foot level where the temperature was -1.5° C and the visibility 
very limited. After 25 minutes of continuously increasing ice accumula- 
tion at this altitude, a slight vibration of the car was noted. About 
5 minutes later this vibration became excessive and caused the instru- 
ment panel to shimmy up to the limits of its shock mounts . The fre- 
quency of the car vibration was estimated as about 5 cycles per second. 
During this period ice thrown from the propellers was heard striking 
the sides of the can which suggested that propeller unbalance may have 
been present. These pieces of ice also caused a rip in the fabric fair- 
ing between the car and envelope. The vibration diminished during the 
following 10 to 15 minutes, then increased again and became very pro- 
nounced about 5 minutes later. A descent was made shortly thereafter 
out of the icing level into above freezing temperatures. The airship 
had been exposed to continuous icing for about 70 minutes. The clear 
ice formed under these conditions is illustrated in figure 3, which 
shows ice formations on the yawmeter mounted below the pilot's window. 

During the descent sheets of ice began sliding off the envelope 
from the bow area and falling back into the propellers and windshield. 
Ice was also noted falling from other areas along the sides of the en- 
velope during the approach to the runway. The surface-ice accretion 
meters were inoperative during this flight, and therefore the envelope 
ice load was not measured. A total ice load of about 4000 pounds was 
estimated by considering airspeed, deck angle, and power settings. This 
included the acc umula ted ice on fin brace wires and protuberances as 
well as on the envelope surface. 

Near the end of the icing period the fin brace wires were observed 
oscillating at an estimated 6-inch amplitude. Ice accretions up to 

1-i inches in diameter had formed on the 3/l6-inch-diameter wires. Some 

flapping motion of the fins was also detected. The control cables had 
also collected some ice, but their operation was not impaired. An in- 
spection upon landing revealed some 30 punctures in the fabric surfaces 
of the two lower fins. These 1- to 6- inch punctures were attributed to 
pieces of ice thrown from the oscillating brace wires or to chunks of 
ice falling from the leading edge of the upper fins or both. 

Other operational problems resulting from the icing included the 
freezing shut of an air damper controlling air to the center ballonet 
and ice on the windshield which completely obscured forward visibility. 
Operation of the damper by either electrical or manual controls was not 
possible. The windshield defogging system was ineffective in preventing 
or removing the ice coating. 
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The car vibration was attributed to propeller unbalance caused by 
uneven throwoff of- ice from the individual blades . A propeller unbalance 
could produce a greater airframe vibration on the ZPG-2 airship than on 
a conventional aircraft because of the propeller mounting configuration. 
Since the propellers are rigidly mounted to the outrigger frames without 
shock mounts, vibrations from unbalance can be easily transmitted through- 
out the airframe. Also, the cantilever design of the low-mass outrigger 
structures allows for marked increases in the amplitude of forced vibra- 
tions under conditions of resonant frequencies. 

Another flight in freezing- rain conditions (flight number 4, table 
i) produced only small ice accretions at temperatures slightly below 
freezing. Intermittent periods of icing existed at the freezing level 
for about 30 minutes. Wires and protuberances collected a slight amount 
of ice. After the airship left the precipitation area, water ran off 
the envelope, which indicated some ice had formed on the top surfaces. 

The ice load was not estimated. 

A vertical sounding was made during this flight to determine the 
depth of the freezing-rain level. Rain at above freezing temperatures 
was falling at the ground at the time of takeoff. Low ceilings (l50 ft) 
and visibilities (l/2 mile) were prevalent in the area of the flight. 
During climb the temperature decreased to freezing at 1700 feet, but at 
3000 feet the temperature was again above freezing. Hence, the layer 
of freezing rain was less than 1500 feet deep. The inversion continued 
to be observed up to the altitude limit of the airship of 5000 feet where 
light rain was falling. 

The only other icing encounter in which freezing drizzle was present 
occurred during flight number 5 (table I). Since the flight was in clouds, 
supercooled cloud droplets may have contributed to the ice formations . 
Effects on the airship were similar to those described for flight number 
3, although total ice accumulations were not as great. In 6 hours of 
flight icing occurred during three separate periods, each approximately 
40 minutes in duration. The encounters were separated by a sufficient 
time lapse (l to 2 hr) at above freezing temperatures to allow most of 
the ice formed during the previous encounter to melt. The ice buildup 
recorded by the icing-rate meters for each icing period is shown in 
figure 4. 

The first period of icing (flight 5, encounter (a), fig. 4) occurred 
at 2400 feet (at -1° C) during climb over South Weymouth in conditions 
described as freezing drizzle . As the maximum buildup of ice was reached 
(l/4 in. recorded by an icing- rate meter), slight car vibration was felt 
which was attributed to propeller icing. Ice also formed on the antennas 
(causing vibrations), on the windshield, on the radorne, and on other com- 
ponents. After about 3/4 hour the shallow freezing layer in the inver- 
sion decreased in thickness to the points where melting began. 
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During the next Itj hours the airship proceeded to the Cape Cod area 

where icing again was encountered. All ice accretions from the previous 
encounter had melted. Vibration again was felt when about l/4 inch of 
ice had accumulated. Excessive antenna whipping was also present. When 
maximum accretion was reached (about l/2 in. in 30 min), there was ex- 
cessive vibration of the car for about 5 minutes . This was accompanied 
by propeller ice throwoff . A climb at this point to 3200 feet produced 
melting at 1° C, which indicated that an inversion was also present in 
this icing area. As the melting continued, ice in thin sheets fell from 
the sides of the envelope and struck the propeller, the outriggers, and 
the windshield. 

The third icing condition was entered upon return to the South 
Weymouth area at 3000 feet. This condition was similar to the first 
encounter of the flight in that only l/4 inch of ice accumulated over 
a 40-minute period. Antenna whipping again appeared excessive. Air- 
speed was reduced to help prevent the wires from snapping. Inspection 
of the airship after landing disclosed ice accretions l/2 to 3/4 inch 
thick on the fin brace wires and control cables, but no ice had remained 
on the envelope surface . 

During a takeoff on the day following this flight, a control problem 
arose which was considered related to the icing encountered on the pre- 
vious day. In an attempt to climb the airship veered sharply to the 
right . Flight was maintained, however, and a landing was effected with 
the controls still fouled . Inspection showed that movement of the lower 
left ruddevator was prevented by the control cable, which had jumped a 
pulley and lodged between the pulley and its guard. The failure was 
assumed to have developed during the rollout phase of the landing from 
the previous flight. A buildup of ice on the affected cable and pulley 
would cause no trouble for normal flight, because the control movements 
are us ually small. However, the large cable travel required on landing 
could cause the enlarged cable to ride up and out of the iced-over pulley 
groove. A similar control-system difficulty was encountered during flight 
number 6 also. When the automatic pilot was disengaged prior to letdown, 
movement of the controls was dangerously restricted. They could be moved 
manually only in s mall increments . Complete control was regained after 
descent into above freezing temperatures . 


Flights in Rime Ice 

The rime icing conditions were found in typical supercooled clouds 
with temperatures ranging from -2° to -10° C and liquid-water contents 
from 0.1 to 0.5 gram per cubic meter. The ice formations were different 
from those produced in freezing drizzle in that the rime accretions were 
confined to the forward edges of small objects only and, unlike glaze- 
ice formations, had less tendency to spread. A comparison of the two 
types of ice formation is shown in the flight photographs of figure 3. 
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The initial flight in rime icing conditions (flight number 1, table 
i) revealed the antenna whipping problem. Of particular concern at the 
time was the possibility of an antenna breaking. Had the starboard an- 
tenna broken near the forward mount, the wire could have tangled with 
the ripcord and whipped back into the propeller, pulling the ripcord with 
it. (The ripcord is an emergency device which when pulled tears an 
opening in the top of the envelope to allow rapid escape of the helium. 
The purpose is to destroy aerostatic lift rapidly in emergency situations 
developing at or near the ground.) 

The second airship, which was flown during the final season of the 
project, encountered only rime icing conditions in operations intended 
to evaluate the ability of the airship to conduct various missions. Test 
objectives were (l) to evaluate extended flights over water (24 hr or 
more) simulating patrol missions in conditions favorable for icing, (2) 
to demonstrate "barrier" station operations for an extended period (10 
days), and (3) to conduct a distance -endurance flight following a cir- 
cuitous route around the northern Atlantic Ocean (8200 miles in 11 days) . 
In the icing conditions encountered during these missions no adverse 
effects on airship operations were noted, and no potentially hazardous 
situations developed. 

This airship was provided with limited icing protection based on the 
flight experiences with the original test airship. Modifications con- 
sidered most desirable and installed on the airship included: 

(1) A shielded control cable system (at ruddevator horn pulleys 
only) 

(2) Flush-mounted antennas 

(3) Electrically heated propellers (external rubber boots) 

(4) Electric ally heated center-ballonet air-inlet dampers 

(5) Increased hot air to the windshield 

(6) Added strength (protective coating) to the fabric covering top 
surfaces of lower fins 

(7) A protective blanket inside the envelope in the plane of the 
propellers 

Flight number 8 (table I) , which was conducted in support of the 10- 
day barrier operation, encountered prolonged periods of rime icing while 
holding in the area. The airship was in clouds and snow during most of 
the 40-hour flight. Synoptic data indicated a situation of mixed clouds 
(supercooled water changing to ice crystals and snow) formed by very 
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cold air over a warm sea surface. Icing occurred in four separate per- 
iods with an estimated cumulative time in icing of about 5 hours. Ice 

accretions (about 1^ in.) were the largest recorded for rime icing en- 
counters. The pilot noted that the ice on small protuberances would 
build up to a certain thickness and break off. Thus, it would seem that 
this condition approached the maximum effect of prolonged periods of 
rime icing. The total ice load was estimated not to have exceeded 2000 
pounds. Some snow in the tail section contributed to the increased 
heaviness . 

In general, the over-all operation of the airship was not seriously 
impaired by the rime icing conditions that were encountered. The limited 
icing protection installed on the second airship appeared adequate for 
those components which were adversely affected by the rime accretions. 


Flights in Snow 

Several flights were conducted in moderate to heavy snow with little 
or no in-flight effects on the airship. Snow loads accumulated while an 
airship is on the mast have long been recognized as a serious problem. 

In flight, however, most of the snow catch blows off at the higher air- 
speeds, except possibly very wet snow. The results of a 2-hour flight 
in light to moderate snow at -2° to -6° C are shown in figure 5(a). A 
pmal 1 catch est imat ed to be about 1000 pounds can be seen distributed 
along the top of the envelope and in the catenary seam along the side. 

On other flights snow was observed only on the top, between and aft of 
the fin area, as shown in figure 5(b). 

During the test flights snow loads created operational problems on 
two occasions. While the airship was on the mast prior to takeoff for 
flight number 9, a heavy snow load collected on the top. Although some 
snow blew off during takeoff, the airship remained very tail-heavy and 
required all ballonet air forward plus the elevator down to effect a 13° 
angle of attack at 2000 feet. The airship was operating at maximum gross 
weight and low power settings to minimize fuel consumption because of the 
endurance requirements of the mission. About 20 minutes of icing was en- 
countered soon after takeoff, which produced a small accumulation (0.1 
in.). The combined effect of all these conditions created a stall situ- 
ation which required descent to a lower altitude. This provided in- 
creased ballonet air capacity, which was used to reduce the high angle 
of attack. With less down elevator also reducing the drag, level flight 
could be maintained without an increase in power. 

Heavy snow produced another problem during a special maneuver con- 
ducted while on a barrier station. An attempt to pick up sea-water bal- 
last was abandoned when altitude could not be maintained at the very low 
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ground speeds required for this operation. The crew concluded that with 
light surface winds (7 to 9 knots) the low airspeed allowed snow to accu- 
mulate, whereas at higher airspeeds the snow would blow off. Another 
consideration was that warmer temperatures near the sea surface produced 
wet snow which adhered to the top surfaces better. 


THEORETICAL ANALYSIS OP AIRSHIP ICING IN FREEZING RAIN 

The flight-test results confirmed the prediction from theory that 
heavy ice loads on the envelope occur only in freezing rain or drizzle. 
Because of the low frequency of occurrence of freezing rain, insufficient 
experience was obtained to provide a reliable basis for estimating max- 
imum icing rates. It appeared desirable, therefore, to resort to calcu- 
lations to determine the relation between icing rate and meteorological 
conditions, in order that climatological data might be used to evaluate 
the operational significance of icing in freezing rain. 


Calculation of Rate of Icing in Freezing Rain 

The rate of ice formation per unit area at a given location on the 
airship envelope depends on the rate of water impingement, the airspeed, 

and the wet-bulb temperature of the air. The temperature dependence re- 

sults from the fact that, as ice is formed, the latent heat must be re- 
moved by convective transfer to the environment. For any combination of 
airspeed, temperature, and location on the airship there exists a crit- 
ical freezing rate at which the release of latent heat is just suffi- 
cient to raise the surface temperature to 32° F. If the water impinge- 
ment rate is less than the critical freezing rate, all the water freezes, 
and the icing rate equals the impingement rate. If the impingement rate 
is greater, most of the excess water runs off in the liquid state. 

In order to determine the total rate of ice formation on the enve- 

lope for a given combination of airspeed, rainfall rate, and wet-bulb 
temperature, it is necessary to calculate the distribution of water im- 
pingement rate and critical freezing rate over the envelope surface . 

The distribution of impingement has been determined as a function of 
rainfall rate and airspeed by assuming that the impingement is the sum 
of two terms: (l) the gravitational impingement that would occur if 

the airship were at rest in still air and (2) the dynamic impingement 
that would result from forward speed alone, neglepting gravity. The 
distribution of the critical freezing rate was calculated from heat- 
and mass-transfer equations based on references 3 and 4. Both sets of 
calculations are described in appendix B (symbols are defined in appen- 
dix A). A comparison of critical freezing rate and impingement rate, 
expressed in pounds per hour per foot along the axis, for a particular 
set of conditions is shown in figure 6(a) . At locations where the 
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critical freezing rate is greater than the impingement rate, all the im- 
pinging water freezes, and the icing rate equals the impingement rate. 
Where the impingement rate exceeds the critical freezing rate, the icing 
rate is greater than the critical freezing rate by the amount of ice 
formation required to release sufficient latent heat to raise the temper- 
ature of the excess water to 32° F. 

For example, in figure 6(a) at a location of 60 feet from the nose 
the impingement rate (33 (lb/hr)/ft) is less than the critical freezing 
rate (44 (lb/hr)/ft), and the icing rate is equal to the impingement rate 
(33 (lb/hr)/ft). At a location 20 feet from the nose, on the other hand, 
the impingement rate (117 (lb/hr )/ft) is greater than the critical freez- 
ing rate (67 (lb/hr)/ft). The excess water (50 (lb/hr)/ft) undergoes a 
temperature rise from 27° to 32° F, which absorbs heat at a rate of 50x5 
or 250 (Btu/hr)/ft, and thus freezes additional ice in the amount of 
250/l44 or 1.7 (lb/hr)/ft. The local icing rate is therefore 67 + 1.7 
or 68.7 ( lb/hr )/ft . 

The total icing rate for the entire envelope is found by integrating 
the local icing rate over the total length. In the coordinates of figure 
6(a), this integration is accomplished by determining the area under both 
the impingement-rate and critical-freezing-rate curves, using whichever 
curve is lower at any location. To this area is added (32 - tj^/lii 

times the area lying below the impingement curve and above the critical- 
freezing-rate curve. 

Figure 6(b) shows how the relation between impingement and critical 
freezing rate changes as the rainfall rate increases at constant temper- 
ature and airspeed. At rainfall rates up to 0.04 inch per hour, for the 
conditions shown in figure 6(b) (28° F, 30 knots), the area of excess im- 
pingement and runoff comprises only a small portion of the envelope at 
the nose. Thus, the icing rate increases almost linearly with the rain- 
fall rate. At a rainfall rate of about 0.047 inch per hour runoff sets 
in along the rear half of the envelope about 250 feet from the nose, and 
at a rainfall of 0.05 inch per hour most of the rear half is running wet. 
Beyond this point the icing rate increases more slowly with increasing 
rainfall. Finally, beyond 0.075 inch per hour the entire envelope is 
subject to runoff, and further increases in rainfall cause only a small 
increase in icing. Curves such as those shown in figure 6 provide a 
means of estimating the total icing rate as a function of rainfall rate 
for a given airspeed and temperature , as shown in figure 7 . 
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Effect of Climatic Factors on Probable Frequency 
and Severity of Icing in Freezing Rain 

The calculated rates of icing presented in figure 7 show that very 
heavy ice loads can be collected in periods of an hour or less in certain 
combinations of temperature and rain intensity. The operational signif- 
icance of these results depends on the frequency of freezing rain and 
the relative frequency of various combinations of temperature and rain- 
fall rate in freezing rain. These are climatic factors, varying with 
geographical location and season. 

The occurrence of freezing rain requires the presence of a low- 
altitude layer of air at temperatures below freezing surmounted by a 
layer of warmer air in which the temperature is above the melting point . 
When a shallow layer of cold air flows over relatively warm water, the 
surface-air temperature rapidly approaches that of the water. Heat- 
transfer calculations described in appendix C show that the existence of 
conditions required for freezing rain at low altitudes over the sea is 
virtually impossible in air that has had a trajectory of approximately 
300 miles or more over water at a temperature of 45° F or higher. In 
the light of this result, an examination of the January sea-surface tem- 
perature distribution (fig. 8) shows that the only potential icing areas 
that need be considered with reference to the aircraft-early-warning 
mission axe the coastal areas of the northeastern United States and the 
adjacent ocean areas . 

Frequency of occurrence of freezing rain in northeastern United 
States . - The estimated average frequency of occurrence of freezing rain - 
in winter in the northeastern United States is shown in figure 9. The 
data shown axe 3-month (December through February) average percentages 
based on two daily observations at 12 stations with an average record 
length of 9 years. These values were obtained from an unpublished list- 
ing of free zing- level and surface -temperature data for freezing precipi- 
tation furnished by the National Weather Records Center. The general 
pattern of freezing-rain frequency shows a gradual increase from south 
to north and a rapid decrease along the coast. This strong coastal ef- 
fect is the result of the warming influence of the ocean in cases in 
which an onshore component exists in the cold-air flow. 

Probable maximum airship icing rates in freezing rain over land area 
in northeastern United States . - Since the icing rate is a known function 
of airspeed, rainfall rate, and free-air temperature, the frequency dis- 
tribution of icing rates for a given airspeed can be calculated from the 
frequency distributions of temperature and rainfall rate. The frequency 
distribution of surface-air temperatures reported during freezing rain 
is shown in figure 10 (data from the same source as fig. 9). Surface- 
air temperatures are probably representative of flight temperatures for 
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aircraft-early-warning airship operation over land at low altitudes dur- 
ing return to or departure from the base. An estimated frequency dis- 
tribution of rainfall rate applicable to freezing rain is shown in fig- 
ure 11. The method of selecting data and the theory underlying their 
use to represent freezing rain are given in appendix D. 

Combining the temperature distribution of figure 10, the rainfall 
distribution of figure 11, and the relation between icing rate and tem- 
perature and rainfall rate shown in figure 7 yields the frequency distri- 
butions of icing rate shown in figure 12 for airspeeds of 30 and 50 knots. 
Thus, if an airship is flying at 30 knots in freezing rain, the chances 
are 1 in 10 that the icing rate exceeds 6300 pounds per hour and 1 in 
100 that it exceeds 14,000 pounds per hour. The temperature and rainfall 
distributions would not be expected to vary widely in the limited geo- 
graphical area of the northeastern United States because the occurrence 
of freezing rain fixes the air-mass temperatures within rather narrow 
limits. Thus, figure 12, which gives the probability of various icing 
rates when freezing rain is known to be occurring, may be regarded as 
applicable generally to land areas near the east coast of the United 
States in winter. The geographical effects within this area are included 
in the frequency of occurrence of freezing rain (fig. 9). 

Estimated conditions off northeastern coast of United States . - As 
implied in figure 9, the probability of encountering freezing rain is 
much less over the ocean off the coast of the northeastern United States 
than it is over adjacent land areas. Since direct observations of 
freezing-rain frequencies at sea are not available, it is desirable to 
estimate the effect of flow over the relatively warm sea surface in 
modifying the air-mass temperature structure in freezing-rain situations. 

An analysis of heat transfer occurring when a layer of cold air 
topped by an inversion flows over a warm water surface is presented in 
appendix C. The distance of travel over the ocean (surface temperature, 
6° C) required to raise the surface-air temperature to 2° C was calcu- 
lated as a function of the initial surface-air temperature and the height 
of the first freezing level. A surface-air temperature of 2° C was 
chosen as representing a condition under which an airship could safely 
descend to below the freezing level (about 700 ft) and shed accumulated 
ice. The relative frequency of various combinations of surface-air tem- 
perature and height of the first freezing level measured during freezing 
rain were combined with the calculated values of distance to 2° C to ob- 
tain the frequency distribution shown in figure 13. This curve shows 
the percentage of cases of freezing rain (observed on land) that would 
be warmed sufficiently to have a surface-air temperature of 2° C or 
higher as a function of distance transported over water at 6° C. Thus, 
one-half of the freezing-rain cases at the coast would be effectively 
eliminated at 64 nautical miles offshore, 90 percent at 145 miles, and 
99 percent at 260 miles . 
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In interpreting these results, it should he noted that synoptic ex- 
perience has shown that most cases of freezing rain on land are not trans- 
ported directly out to sea. Thus, figure 13 overestimates the hazard, 
since it is based on the assumption that all freezing- rain situations 
move out to sea. 


CONCLUDING REMARKS 

In test flights conducted by the Navy to evaluate the effects of ice 
accretions on the ZPG-2 airship a number of typical icing conditions were 
encountered. While the airship was exposed to these conditions, ice 
formed on various components but presented no serious hazard to the oper- 
ation of the airship, although sane potentially unsafe situations devel- 
oped. The more common rime ice formed only on the forward edges of 
small protuberances and wires without causing appreciable concern. Ice 
accretions of the glaze type which occurred in what was described as 
freezing drizzle adversely affected various components to a somewhat 
greater extent. Protective designs for some of the affected components 
appeared desirable to assure safe operations. 

A study of the ice accretion problem for airships indicated that 
in addition to component icing a problem could develop from an excessive 
ice load on the large surface area along the top of the envelope which 
would occur only in conditions of freezing rain. This potentially haz- 
ardous situation was not adequately investigated in the flight tests be- 
cause of the infrequent occurrence of these conditions. However, the 
tests verified that ice can form on the top of the envelope when the 
drop sizes approach those characteristic of freezing rain. During these 
flights the ice formations did not produce a sufficient increase in air- 
ship heaviness to become an operational problem. 

The calculations presented in this report show that very heavy ice 
loads can be collected under more extreme conditions of rain intensity 
and temperature than those encountered during the flight program. The 
operational significance of these hazardous conditions depends on the 
frequency of their occurrence. In a search for all types of inclement 
weather the flight test program demonstrated the low probability of these 
situations occurring over land areas along the northeastern coast of the 
United States. The probability is further evaluated in this report using 
several years of meteorological data from land areas . 

The calculated icing probability curves and the calculated effect 
of warming over the sea are somewhat speculative. Actual numerical val- 
ues are therefore uncertain, but the following qualitative conclusions 
of operational significance may be stated with considerable confidence: 
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(1) It is possible to encounter very heavy rates of icing under 
conditions of high rainfall rate and low temperature. 

(2) Statistics on rainfall rate and temperature in freezing rain 
indicate that combinations giving rise to hazardous icing rates may be 
expected in from 1 to 5 percent of freezing-rain occurrences over land 
areas in the northeastern United States. Since freezing rain occurs in 
this area in only about three to five storms per year, the probability 
of severe icing is very small. 

(3) Because of the warming effect of the sea surface, the probability 
of encountering hazardous icing in freezing rain decreases rapidly with 
distance offshore, becoming negligible at 200 to 300 miles. 

(4) The effect of freezing rain in aircraft-early-warning operations 
is likely to be confined to arrivals and departures, when the airship is 
over land or a short distance offshore. Hazardous icing during these 
phases of operation can probably be avoided with the aid of weather fore- 
casts, since the synoptic conditions required for freezing rain can be 
forecast satisfactorily for short periods. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, December 17, 1957 
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APPENDIX A 
SYMBOLS 

frontal area of airship, sq ft 

horizontal projected area of airship, sq ft 

specific heat of air at constant pressure, Btu/(lb mass)(°F), 
kcal/(kg) (°C) in appendix C 

diameter of envelope at distance z from nose, ft 

collection efficiency, dimensionless 

critical freezing rate per unit distance along axis, lb mass/ 
(hr) (ft) 

critical freezing rate per unit area, lb mass/(hr)(sq ft) 

heat-transfer coefficient, Btu/(hr)(sq ft) (°F) , kcal/(hectare) 
(hr) (°C) in appendix C 

convection heat-transfer coefficient, kcal/ (hectare) (hr) (°C) 
(appendix C) 

radiation heat-transfer coefficient, kcal/ (hectare) (hr) (°C) 
(appendix C) 

specific enthalpy of ice at 32° F, Btu/lb mass 

specific enthalpy of liquid water at kinetic wet-bulb tempera- 
ture, Btu/lb mass 

specific enthalpy of liquid water at wet-bulb temperature, 

Btu/lb mass 

specific enthalpy of water vapor at kinetic wet-bulb temperature, 
Btu/lb mass 

specific enthalpy of water vapor at 32° F, Btu/lb mass 
mechanical equivalent of heat, ft-lb/Btu 
thermal conductivity of air, Btu/(br) (ft) (°F) 
total rain impingement rate, lb mass/hr 
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M d 

M g 

m d 

tn g 

Nu 

P 

Pr 

p kw 

P 0 

P 32 


%s,e 


d as,fr 


a s s , c 


ss , e 


*ss,r 


R 

Re 

r 

st c 

St, 


cond 


dynamic impingement, lb mass/hr 

gravitational impingement, lb mass/hr 

dynamic impingement per unit length, lb mass/ (hr) (ft) 

gravitational impingement per unit length, lb mass/ (hr) (ft) 

Nusselt number, dimensionless 

atmospheric pressure, mb 

Prandtl number, dimensionless 

saturated vapor pressure over liquid water at t^, mb 
free -stream vapor pressure, mb 
saturation vapor pressure at 32° F, mb 
heat of vaporization of water, kcal/kg (appendix C) 
increase in enthalpy of cold-air layer, kcal/(hectare) (hr) 
heat loss from airship surface by convection, Btu/(hr)(sq ft) 
heat loss from airship surface by evaporation, Btu/(hr)(sq ft) 
heat liberated at airship surface by freezing, Btu/(hr)(sq ft) 
heat loss from sea surface by convection, kcal/ (hectare) (hr) 
heat loss from sea surface by evaporation, kcal/ (hectare) (hr) 
heat loss from sea surface by radiation, kcal/ (hectare) (hr) 
r ainf all r at e , in . /hr 

Reynolds number based on surface distance s, dimensionless 

temperature recovery factor, dimensionless 

Stanton number for convection, dimensionless 

Stanton number for radiation, dimensionless 

distance along airship surface from nose, ft 

temperature at condensation level, °C 
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t; 


cond 


T- 


L o,cl 


T 


°,y 


ss 


ss 


'ki 


D kw 


0 


U 


U, 


0 


v 

w 

x l 

x t 

*v 


■S', sat 

y 

z 

^cond 

"Tr 


temperature at condensation level, °K 

temperature on dry adiabatic line through T c , °C 

surface-air temperature, °C 

surface -air temperature at coast line, °C 

surface -air temperature after trajectory of length y, °C 

sea- surface temperature, °C 

sea-surface temperature, °K 

temperature in cloud above condensation level, °C 

kinetic iced-bulb temperature, °F 

kinetic wet-bulb temperature, °F 

free-stream temperature, °F 

local airspeed outside boundary layer, ft/sec 

free-stream velocity, ft/sec 

airspeed or wind speed, knots * 

raindrop falling speed, knots 

liquid-water content in rain, lb mass/(sq ft) (nautical mile) 
liquid-water mixing ratio, g water/g air 
total water mixing ratio, g water/g air 
water-vapor mixing ratio, g water/g air 

saturation water-vapor mixing ratio over sea water at T g 

distance over ocean, nautical miles 

height above sea surface, hectometers 

height of condensation level, hectometers 

height of freezing level in inversion, hectometers 
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a 


e 

X 

M- 

P 

p , 

P 

a 

T 

§ 


height of inversion base, hectometers 
pressure altitude, ft 

distance from nose, measured along axis, ft 

axial distance to limit of dynamic impingement, ft 

ratio of dynamic to gravitational impingement, dimensionless 

depression of kinetic wet-bulb temperature below freezing, °F 

wetted-area factor, dimensionless 

viscosity of air, lb mass/ (ft) (sec) 

density of air, lb mass/cu ft, kg/ (hectometer)^ in appendix C 
density of air at sea level, kg/ (hectometer) 

average air density from sea level to height of inversion base, 

kg/ (hectometer) ^ 

Stefan -Boltzmann constant, kcal/ (hr) (hectare) (°K) 4 
t ime , hr 

thermodynamic heat-transfer function defined in eq,. (Bll), 
lb mass/(hr)0’2(ft)l*8(nautical mile)0.8 
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APPENDIX B 


CALCULATION OF AIRSHIP ICING RATE IN FREEZING RAIN 

The calculations presented herein apply to the "theoretical envelope" 
only and do not include icing of the car, fins, brace cables, and so 
forth. 


Impingement of Rain on Airship Envelope 

In estimating the rate and distribution of rain impingement on an 
airship it is necessary to consider two processes, gravitational settling 
and dynamic interception. These two processes are not independent, be- 
cause gravitational settling is modified by the local airflow about the 
airship, and the dynamic impingement is influenced by the effect of gravity 
on the raindrop trajectories. Because of the difficulties involved in an 
analytical treatment of combined gravitational and inertial impingement, 
and since trajectory calculations neglecting gravity were at hand, the 
two effects were evaluated independently, with the assumption that the 
total impingement is equal to two terms: (l) the gravitational impinge- 

ment Mg that would occur if the airship were at rest in still air, and 
(2) the dynamic impingement M d that would occur if the airship were 
moving through a field of raindrops suspended in air in the absence of 
gravity . 

The gravitational impingement is a function of the rainfall rate and 
the projected area. Since 1 inch of rain equals 5.2 pounds per square 
foot, the total gravitational impingement is given by 

M g = 5.2 RA h (Bl) 

The gravitational impingement per foot of length measured along the axis 
of the envelope is 


m g = 5.2 RD Z (B2) 

The distribution of gravitational impingement along the length of the 
airship is represented by m g /R (lb/(ft)(in. rain)), as shown in figure 14. 

The dynamic impingement depends primarily on drop size, liquid-water 
content, airspeed, and shape and size of the airship. The equation for 
dynamic impingement is 


M d = wVA^E 


(B3) 
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where the collection efficiency E is the fraction of drops in the path 
of the body that strike the surface. Calculated values of E are avail- 
able for spheres and prolate ellipsoids of revolution of fineness ratios 
5 and 10 (refs. 5 to 7). Since the portion of the ZPG airship envelope 
forward of the maximum- thickness point is very nearly an ellipsoid of 
fineness ratio 3.25, values of collection efficiency were obtained by 
interpolation. The results are shown in figure 15 in terms of drop size 
and airspeed. The values shown were determined for a 3000-foot altitude 
and -2° C and are applicable approximately from sea level to 5000 feet 
and from 0° to -10° C. 


The liquid-water content may be eliminated from equation (B3) by 
introducing the rainfall rate and falling speed: 

w = 5.2 R/v 


V 


M d = 5.2 RAp ^ E 


The ratio of the dynamic to the gravitational term is 

Af v 


M d 

— = a = -r— — E 


M, 


g 


A h v 


and the total impingement is 

M = 4* Mg — 5*2 (l 4* cl) 


For the ZPG airship. 


M = 1.02X10 5 R(1 + a) 


(B4) 


(B5) 


(B6) 


(B6a) 


The falling speed of raindrops at 3000 feet is given as a function of 
drop diameter in figure 16 (based on refs. 8 and 9). Values of a for 
the ZPG airship, calculated from the data of figures 15 and 16, are pre- 
sented in figure 17 as a function of drop diameter and airspeed. 


Th^ distribution of dynamic impingement as a function of the distance 
z measured along the axis from the nose was determined from the data of 
references 5 to 7. It was found that the local dynamic impingement rate 
in pounds per hour per axial foot is a linear function of z as given by 


m d 




(B7) 


where z ffl is the axial distance from the nose to the limit of dynamic 
impingement. Values of z m obtained by interpolation from references 5 
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to 7 are shown in figure 18 as a function of drop diameter and airspeed. 
If the local dynamic impingement rate at the nose is divided by the rain- 
fall rate, the result is 


1 V, q ~ 

— - - 2.04X10° — 
z m R z m 


(B8) 


which represents the impingement on the nose in pounds per axial foot per 
inch of rain. This quantity is a function of airspeed and drop diameter 
and is shown in figure 19. 

Figures 15, 17, 18, and 19 define impingement as a function of drop 
size and airspeed. For a given rainfall rate, the collection efficiency 
increases and the liquid-water content decreases with increasing drop 
size. As a result, the variation with drop size of total impingement is 
relatively small for drops greater than 0.6 millimeter in diameter. As 
shown in figure 17, the quantity (1 + a) changes by only 9 percent at 30 
knots and 11.5 percent at 50 knots for a change in drop diameter from 0.6 
to 1.8 millimeters. For this reason, no serious errors are introduced by 
using a statistical average relation between rainfall rate and drop diam- 
eter to express the impingement in terms of rainfall rate and airspeed. 
This relation, shown in figure 20, was derived from an empirical equation 
relating rainfall rate and liquid-water content from reference 10. 

The elimination of drop size as an independent variable gives the 
following impingement results for the ZPG airship: 

(1) Total impingement (lb/hr) as a function of airspeed and rainfall 
rate (fig. 21) (The gravitational impingement alone is shown by 
the curve for zero airspeed.) 

(2) The local dynamic impingement rate at the nose (lb/ (hr) (axial 
ft)) as a function of airspeed and rainfall rate (fig. 22) 

(3) The limit of the dynamic impingement zone (ft from nose along 
axis) as a function of airspeed and rainfall rate (fig. 23) 

The distribution of total impingement along the axis may be found 
by adding the local gravitational impingement, from figure 14, to the 
local dynamic impingement, obtained from equation (B7) using z from 
figure 23 and 2M^/z m from figure 22. 

A convenient rule for estimating the total impingement rate which 
gives results within 5 percent of the values shown in figure 21 is as 
follows: The total impingement in thousands of pounds per hour is equal 

to the rainfall rate in hundredths of an inch per hour increased by 1 
percent for each knot of airspeed in excess of 10 knots. 
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Effect of Heat Transfer on Icing Rates 


In conditions of light freezing rain at temperatures considerably 
below freezing, all the water impinging on the airship freezes. Thus, 
the rate and distribution of icing are the same as the rate and distribu- 
tion of impingement. As the water freezes, the liberated heat of fusion 
causes the surface temperature to rise above the air temperature until 
the rate of heat loss from the surface by convection and evaporation 
balances the rate of liberation of latent heat. At higher rates of water 
impingement the surface temperature may reach 32° F. Additional increases 
in impingement cannot increase the surface temperature above 32° p and 
therefore do not increase the evaporation and convection heat losses. 

The increase in ice formation due to an increase in impingement under 
these conditions is only enough to provide sufficient latent heat to 
raise the temperature of the additional impinging water to 32° F. The 
remainder of the additional water runs off in the liquid state. For any 
particular combination of airspeed, altitude, temperature, humidity, and 
position on the airship, there exists a critical rate of water impinge- 
ment which is just sufficient to maintain the surface temperature at 32° 

F. The rate of ice formation per unit area under these conditions is the 
critical freezing rate f cr . 


Under these conditions, the rate of heat loss from the surface by 
convection q oc . . and evaporation q„„ is balanced by the rate of 

ao y G db ^ c 

production of heat by freezing q ag £ r . The rate of heat loss by radia- 
tion is much smaller than the convection term and is therefore neglected. 
Also neglected is transfer of heat into or out of the underlying surface. 
The convection term is: 


1 


as,c 




(B9) 


The evaporation term q is written to include the entire heat loss 

as y g 

due to the change of enthalpy of the water that evaporates; the initial 
temperature of the rain drops is assumed to be the wet-bulb temperature: 


q m JL 

^■as.e c-n 


0.622 (p^2 " Pq) 


( u n\ 

(h __ - h, - ^7 ) (BIO) 

\v,32 Z, w 2J / ' ' 


The rate of gain of heat due to the change of enthalpy of the water that 
freezes is 


1 as,fr 


(. 


h + 
Z,w 


u; 


_o 

2J 


" h i,32) 


(Bll) 
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Use of the kinetic wet-bulb temperature t^ w as a reference tempera- 
ture is desirable because it makes the results applicable to unsaturated 
as well as saturated ambient conditions. The kinetic wet-bulb temperature 
is defined as the equilibrium temperature of a wet-bulb thermometer in a 
moving airstream when the wet bulb is covered with liquid water. When the 
airship is flying in clouds, the reading of a dry-bulb thermometer shielded 
from direct droplet impingement may be used as tj^. When in clear air 
at temperatures below freezing, the wet-bulb thermometer yields the kinetic 
iced-bulb temperature t^^, because the bulb is coated with ice instead of 
liquid water. The corrections shown in figure 24 may be used to obtain 
t^ w from t^ and the relative humidity. 

Equations (BIO) and (Bll) may be simplified by using the following 
approximation : 

h Z, w + 2J ” *H,kw 


The kinetic temperature rise of the water drops is about half the wet- 
bulb kinetic temperature rise, and both are less than 1° F at speeds up 
to 60 knots. Thus, the error in the approximation is small compared with 
the enthalpy changes in evaporation (eq. (BIO)) and freezing (eq. (Bll)). 
With this simplification, the terms are substituted in the heat-balance 
equation q as ^ c + q as ^ e = q as ^ fr , and the result is solved for f cr : 


f cr = 




kw 


H 

nr 


1,32. 




32 - t„ - 


rU, 


0 


0 ' 2Jc" 


0.622 


( p 52 ~ Pq) 


(h 


v,32 


- h, . ) 
Z,kw' 




(B12) 


The psychrometric equation defining the kinetic wet-bulb temperature 
with respect to liquid water is 


^kw - ^0 + 


2 

rU 0 0.622 (Pkw ~ p p) 


2Jc 


P 


(■* 1 v,kw " •* 1 Z,kw) (BPB) 


Combining equations (B12) and (B13) and rearranging the vapor-pressure - 
enthalpy products give 


*cr 


H 


^ 1 Z,kw _ ^ 1 i,32 


32 - t 


kw 


0.622 

Pc 


[(p 


32 * p kw 


^ ^v,32 ” ^Z,kw^ 


■fc 


3Z)}\ 

(B14) 


(PO " Pkw^ (\-,kw~^v, 



NACA TN 4220 


25 


Since (p Q - p kw ) is less than (p 32 - pj and (h^ - \ ^?) is less 
than 0.5 percent of ( 1^33 - h^^) for 32>t kw >20, the product 
(p Q " p kw) kw " h v,32^ can be ^glected in equation (B14); hence. 


f cr " 


JL 


H, kw 


h,- 


i ,32 


32 


^kw + 


0.622 

Pc„ 


(P 


32 


Pkw^ ^v,32 



This relation gives the critical freezing rate in terms of the local hea't- 
transfer coefficient and the kinetic wet-bulb temperature. Because of the 
large size of an airship, the boundary layer is' turbulent at ordinary 
flight speeds except for a very small area at the nose. The heat-transfer 
coefficient is given by the following equation, which gives good agreement 
with experimental data for blunt bodies of revolution in turbulent flow 
(ref. ll) : 

Nu = ~ = 0.0296 Pr 1 / 3 Re 0 ' 8 (B16) 

K 


which gives 

H = 0.0296 k Pr 1 / 3 (p/p ) 0,8 U ^’ 8 s " 0,2 (B17) 

With Ug * 8 = (U s /U 0 )°‘ 8 (1.689 V )°* 8 = 1.521 (Ug/Uo) 0 '^ 0 ’ 8 , substitution 
in equation (B15) gives 

f cr =$ CU S /U 0 )°- 8 V°- 8 /s °- 2 (B18) 

where 


0.0450 kPr 1 / 3 (p/p) 0,8 . 0.622 

* = h Z,kw - *i,32 L ‘ ^ ^7 


( p ; 


32 


Pkw^N-,32 “ h Z,kw^j 

(B19) 


is a function of altitude and kinetic wet-bulb temperature. Values of $ 
were calculated for a range of pressure altitudes from 1000 to 5000 feet 
and t^ from 32° to 16° F. The quantity kPr 1 / 3 (p/p ) 0-8 was evalu- 
ated at t kw . Values of U s /Uq were calculated from surface -velocity 
data for an ellipsoidal body from reference 11 . 

It was found that $ could be represented within ±Q.5 percent over 
this range of pressure altitude and temperature (1000 to 5000 ft and 16° 
to 32° F) by the following approximation 
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<t> = (0.0080 0 + 0.000025 0 2 )(l - 0.0000143 Zp) (B20) 

where 0 = 32 - is the depression below freezing of the kinetic wet- 

bulb temperature. 

Over a limited range of temperature from 32° F to about 20° F, $ is 
almost a linear function of 0 . For convenience and simplicity in the 
application of critical-icing-rate calculations, the following linear 
approximation was used in the calculations presented in this report: 

* = 0.0082 0 (1 - 0.0000143 Zp) (B2l) 

For altitudes between 1000 and 5000 feet, the error arising from using 
equation (B2l) instead of equation (B19) is approximately 1 percent at 
0 = 4, 0 at 0=8, -1.3 percent at 0 = 12, and -2.8 percent at 0 = 16. 
Combining equations (B2l) and (B18) gives the final expression for the 
critical freezing rate per unit area: 

f cr = °- 0082 9 (1 - 0.0000143 Zp) ( U s /U 0 ) 0,8 V 0,8 /s°- 2 

The critical freezing rate per unit distance along the axis of the 
airship is given by 


F d) ^ 

cr cr z 2 dz 


where X is the wetted-area factor. In the region next to the nose, 
where the envelope is wetted all the way around, X = 2; aft of the dynamic 
impingement area, where the top half only is wetted, X = 1. The transi- 
tion from X = 2 to X = 1 occurs over a small range of z, just forward 
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APPENDIX C 


CALCUIATION OF OFFSHORE HEATING OF COLD AIR 

It is required to calculate the time rate of increase of temperature 
T q at the base of a cold-air layer of depth Z ib flowing over a rela- 
tively warm water surface. The bottom of the layer (Z = 0) is taken as 
outside the boundary layer at about 8 meters above the sea surface . 
"Surface" temperature and wind speed refer to this height. It is assumed 
that the cold-air layer is thoroughly mixed at all times and thus has 
constant total water mixing ratio x^ throughout, a dry adiabatic lapse 
rate up to the condensation level ^’cond-’ an< ^- a saturated adiabatic lapse 
rate from Z cond to Z ib . If the dry adiabatic line is extended from 



z cond Z ib> ttie "temperature represented by a point on this line at 

height Z is the dry adiabatic temperature T d . The corresponding tem- 
perature T w on the saturated adiabatic line is greater than T d be- 
cause of the latent heat released by condensation of liquid water: 

Qv 

^w = ^d + c X Z 
P 


(Cl) 
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Since this analysis applies to temperatures near and below freezing, the 
total water mixing ratio x t = xj + x v is always less than about 0.005. 
Thus, the specific heat and density of the mixture of air, vapor, and 
liquid water may be taken as approximately equal to those of air at the 
same pressure and temperature. 


When heat and water vapor are introduced and mixed throughout the 
layer,, the time rate of increase in enthalpy per unit horizontal area is 



From equation ( Cl) (for the interval Z con( ^ to Z^) 

dT w _ ^d + ^v 
dx ~ dx Cp dx 



Since the dry adiabatic lapse rate is constant, dT^/dx is independent 
of Z and equal to dT Q /dx as long as x^ = 0 at Z = 0. Also, the 
total mixing ratio x^ is assumed independent of Z; therefore, 



(C4) 

The rate of enthalpy loss from the sea surface per unit area 
of a convection term^ 

is the sum 

^■SS;C = ®c^ss ” *o) 

(C5) 

an evaporation term. 


^ss^e “ ( H c/ C p) ^v^ x v ; sat x t^ 

(C6) 

and a radiation term, 


^-ss^r = a ( T ss ~ T cond^ = ^r^ss " T cond^ 

(C7) 
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Equation (C7) is based on .the assumption that the sea: surface and the 
cloud base radiate essentially as black bodies at temperatures near 
275° K. 


: The heat-balance equation is 


q = q + q + q 

a ss.,c ss,e ss,r 


Thus, 

eLT fly 

^ib pc pdT + ^ib pi ^v dx = ^c^ss ” ^o^ + ^v^^^sat “ + H r ('T sg ” ^cond^ 


(C8) 


The rate of increase of water content in the cold-air layer equals the 
rate of evaporation from the surface; therefore. 


Z 


_ ^t 
ib p dx = 


H 


v,sat 



(C9) 


Subtracting equation (C9) from equation (C8), dividing by c_p , and in- 
troducing the wind speed V = dy/dx give * ° 


Z. 

ib p Q dx 


In 


Vo V 


( t sb- To) If + 


H T 


dx ‘ c„p 0 V 


(^ss “ ^cond) 


dy 

dx 


(CIO) 


Introducing the Stanton number St = H/l8.53CppV (one knot = 18.53 
hectometers/hr) and solving for dT 0 /dy give 

dT 0 Prt 18.53 •- 

dy = z il3 C«c(T„ - T 0 ) + St r (T ss - T con(i )j (Cll) 

Use of this equation for numerical calculations requires values of the 
Stanton numbers for convection and radiation. The Stanton number for 
convection heat transfer was obtained from analyses of the turbulent 
transfer of water vapor at the sea surface (refs. 12 and 13). The prin- 
ciple of the similarity of heat and mass transfer in turbulent flow al- 
lows the calculation of the Stanton number from the evaporation coeffi- 
cients. . For a rough surface (V > 13 knots), the theoretical value of 
the Stanton number is 0.0031 based on a roughness length of 0.6 centi- 
meter and a height of 8 meters for the measurement of surface temperature 
and wind. For a smooth. surface, the theoretical value varies slightly 
with the wind speed but is close to an average of 0.0010 for winds from 
4 to 13 knots. A study of average evaporation over world ocean areas 
(ref. 12) yields an empirical average evaporation coefficient corresponding 
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to a Stanton number of 0.0020, about midway between the theoretical val- 
ues for rough and smooth surfaces. Evaporation coefficients calculated 
from vapor-pressure-gradient data given in reference 13 correspond to 
St c values from 0.0004 to 0.0019 at winds less than 12 knots and St c 
values from 0.0023 to 0.0034 at winds over 16 knots. The values of St c 
used herein are the theoretical values. For V < 13 knots. 


For V > 13 knots. 


St c = 0.0010 


St c = 0.0031 

The Stanton number for radiation was calculated from equation (C7) 
using an average value of T* Qn(i of about 271° K and T^ g = 279° K. 
When V is in knots, "' L 

St r = 0.0072/V 


In order to apply equation (Oil) to calculate the probable effect 
of ,a trajectory over water on a given freezing-rain air -mass structure, 
it is necessary to know the temperature as a function of height for the 
initial conditions when the air crosses the coast. The summary of tem- 
perature and freezing-level data in freezing rain referred to previously 
in the discussion of figure 9 includes a listing of the surface tempera- 
ture T Q in °G and the height above the ground of the first freezing 
level Z^. in hectometers. These two quantities were used to define an 

ideal temperature-height curve as follows. A lapse rate of 0.72° C per 
hectometer was assumed from the surface to a height of l/2 Zf** . A 


linear increase in temperature was assumed from this point to the freez- 
ing level at Zfy. Comparison with actual soundings shown in figure 25 


shows that the assumptions are realistic. This relation is expressed in 
the following equation, in which Z fr and Z^ are in hectometers and 
Tq and T 0 ^ cZ in °C: 


Zib = 



(l + — T °^ CZ ) 

V °* 72 Zft- - T o,c iJ 


(C12) 


The surface relative humidity in cold polar air at sea is generally from 
70 to 80 percent (ref. 14), corresponding to a condensation level of 
about 400 meters. Thus, a constant value of 4° C was assumed for 
T 0 - Tccmd* A constant value of 6° C was chosen for T ss , as represent- 
ing January conditions offshore at the latitude of Boston (ref. 15). 

The value of Z^ from equation (C12) and the values T ss = 6 and 
T cond = - 4 were substituted in equation (Cll), and the result was 
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integrated with respect to T 0 , in order to determine the distance off- 
shore y (nautical miles) at which T attains a partic ul a r value T Q v 



(C13) 


where 

a = (0.72 Zfj. - T 0 ^ cZ )(6 St c + 10 St r ) 


h = (0.72 Zfj. - To,cl)(®^c t St r ) 


c — 0.72 Zfy.- 2T 0 ^ C £ 

In figure 26 is shown the calculated distance offshore y at which the 
surface -air temperature reaches 2° C for winds of 5, 10, and 20 knots. 

The distance y is given as a function of the initial surface-air tem- 
perature at the coast T 0}C l and the height of the freezing level Z^ 
assuming the sea-surface temperature is 6° C. The numbers between the 
curves in figure 26 are numbers of observations of temperature nnri freez- 
ing level (measured over land) in the intervals indicated. The variation 
of offshore heating with wind speed is a result of the variation of St c 
and St r . The convection heat transfer per mile increases by a factor of 
about 3 at the critical wind speed (about 13 knots) when the surface 
changes from smooth to rough. The heat loss from the sea per hour by 
radiation is independent of windj therefore, the radiation transfer per 
mile varies inversely with the wind speed. Thus, the total heating per 
mile approaches a maximum at very low wind speeds with radiation pre- 
dominant. At winds just below 13 knots the heating per mile is a mini- 
mum with convection accounting for about two-thirds of the total. At 
winds over 13 knots, convection heating per mile is constant at a rela- 
tively high value, while the radiation effect slowly approaches zero with 
increasing wind. 

In the absence of actual data on wind frequencies, it was assumed 
that the median wind speed is about 13 knots, with half the cases rep- 
resented by the results calculated for 20 knots and the other half 
equally divided between 5 and 10 knots. Thus, the frequencies in figure 
26 were weighted in the ratio 1:1:2 for winds of 5, 10, and 20 knots and 
combined in figure 13 to show the percentage of cases of freezing rain 
that would be eliminated by heating to 2° C as a function of length of 
trajectory over water at 6° C. 
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APPENDIX D 


DERIVATION OF RAINFALL FREQUENCY DISTRIBUTION 
APPLICABLE TO FREEZING RAIN 

The frequency distribution of rainfall rate in freezing rain is 
difficult to establish because statistics on rainfall rate as ordinarily 
presented do not include simultaneous data on temperature and state 
(solid or liquid) of the water particles. From theoretical considera- 
tions it is known that, for a given amount of lifting, the intensity of 
precipitation depends on the moisture content, and thus also on the tem- 
perature of the air that is being iifted. Rainfall rates therefore in- 
crease on the average with increasing temperature of the air mass in 
which the rain is formed. In the case of freezing rain, the temperature 
of the base of the warm air mass in which the rain is formed, is usually 
just a few degrees above freezing regardless of the cold-air-mass tem- 
perature. Thus, a reasonable estimate of the rainfall intensity distri- 
bution in freezing rain may be obtained from the distribution of ordinary 
rain intensity in the temperature range from about 32° to 45° F, and this 
rainfall distribution may be assumed to.be independent of the temperature 
in the cold-air layer. " 

The climatological data published by the U.S. Weather Bureau (refs. 

16 and 17) contain hourly precipitation totals and daily maximum and 
minimum temperatures. The desired rain intensity distribution was ob- 
tained by selecting days on which precipitation was general in New England 
and tabulating the hourly precipitation amounts recorded at stations hav- 
ing minimum temperatures of 32° F or higher and maximum temperatures of 
45° F or lower. The resulting rainfall frequency distribution is shown 
in figure 11. 
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TABLE I. - DATA FROM ZPG-2 AIRSHIP ICING FLIGHTS 
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Figure 1. - Side view of ZPG airship. (Obtained from U.S. Navy.) 
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Figure 2. - Measured ice accumulation and associated effects occurring during test flight n um ber 
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C-46687 


Airflow 


(a) Clear ice formed during flight number 3 in freezing drizzle. 

Figure 3. - In-flight photographs of exposed yawmeter showing ice accretions 
formed under free zing- drizzle and rime- ice conditions. (Obtained from U.S. 
Navy. ) 
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(b) Rime-ice formed during flight number 8. 

Figure 3. - Concluded. In-flight photographs of exposed yawmeter showing ice 
accretions formed under freezing- drizzle and rime-ice conditions. (Obtained 
from U.S. Navy.) 



Ice accumulation, in 


NACA TN 4220 


39 





0 8 16 24 32 40 48 56 64 72 

Lapsed time in icing, min 


Figure 4. - Measured ice accumulation and associated effects occurring during test flight number 5 
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(b) Snow in aft area between top fins following two 40-hour flights. 

Figure 5. - Concluded. Location of snow catch observed on airship envelope. (Obtained from U.S. Navy.) 
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Figure 6. - Comparison of critical freezing and impingement rates. Pressure altitude, 3000 feet. 
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Figure 7. - Icing rate as function of rainfall rate for various temperatures. Pressure altitude, 3000 feet. 
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Figure 7. - Concluded,. Icing rate as function of rainfall rate for various temperatures. Pressure altitude, 3000 feet. 
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Figure 8. - Average January sea-surface temperatures (°F) near coasts of North America. (Based on data 
from ref . 15 . ) 
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Figure 9. - Percentage of observations in northeastern United States in which 
freezing rain was reported. Averages for December , January , and February 
based on two daily observations for periods of 7 to 14 years. Limits of 90 
percent confidence band shown in parentheses. 
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Figure 10. - Frequency distribution of surface -air temperature in freezing rain at 
selected stations in northeastern United States. 
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Figure 12 . - Probability of exceeding various icing rates in northeastern United States 
during winter when flying in freezing rain. 
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Figure 13. - Calculated heating effect of sea surface at 6° C on freezing-rain air- 
mass temperature. 
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Figure 15. - Collection efficiency of airship in rain. Values calculated for ellipsoid of fineness ratio 
3.25 with maximum diameter of 75.4 feet; pressure altitude, 3000 feet; temperature, -2° C; gravity 
neglected. 
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Figure 16. - Falling speed of raindrops at 3000 feet 
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Figure 17. - Ratio of dynamic to gravitational impingement. Pressure altitude, 3000 feet; temperature, -2° 
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Figure 18. - Limit of dynamic impingement as function of raindrop diameter. 
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Figure 22. - Local dynamic impingement at nose 
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Figure 23. - Limit of dynamic impingement as function of rainfall rate. 
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Figure 25. - Examples of radiosonde observations during freezing rain. 
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Surface-air temperature at coast; °C 


(a) Wind speed, 5 knots. 

Figure 26. - Distance of flow over water at 6° C to attain surface-air temperature of 
2° C as function of initial surface-air temperature and height of first freezing level. 
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(b) Wind speed, 10 knots. 

Figure 26. - Continued. Distance of flow over water at 6° C to attain surface-air 
temperature of 2° C as function of initial surface-air temperature and height of 
first freezing level. 
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(c) Wind speed, 20 knots 

Figure 26. - Concluded. Distance of flow over water at 6° C to attain surface-air 
temperature of 2° C as function of initial surface-air temperature and height oi 
first freezing level. 


NACA - Langley Field, Va. 


